Cooperative scattering measurement of coherence in a spatially modulated Bose gas 
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Correlations of a Bose gas released from an optical lattice are measured using superradiant scat- 
tering. Conditions are chosen so that after initial incident light pumping at the Bragg angle for 
diffraction, due to matter wave amplification and mode competition, superradiant scattering into 
the Bragg diffracted mode is preponderant. A temporal analysis of the superradiant scattering gain 
reveals periodical oscillations and damping due to the initial lack of coherence between lattice sites. 
Such damping is used for characterizing first order spatial correlations in our system with a precision 
of one lattice period. 
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The realization of Bose-Einstein condensation has ren- 
dered possible studies of completely coherent interactions 
between light and matter waves. The observation of su- 
perradiant scattering with a Bose-Einstein condensate 
(BEG) in an exemplary experiment by Ketterle et al. 
[l| lead later on to the demonstration of amplification of 
matter maves In such cooperative scattering pro- 

cesses, light scattering and subsequent atomic recoil are 
deeply modified and tend to be coherent and directional. 
A method based on Raman superradiant scattering was 
recently used to analyse the coherence properties of a 
BEC [5|. Exploiting the strong coherence dependence of 
the light and matter waves mixing occurring in superra- 
diance [sl, Q , it was able to spatially discriminate between 
normal and superfluid phases of the condensate. 

No further studies have been carried out yet to esti- 
mate whether or not superradiance can be a good candi- 
date for spatially characterizing coherence in more com- 
plicated or correlated systems formed after loading Bose 
gases into optical lattices. In principle, coherence proper- 
ties can be accessed directly through interferences imaj 
ing. They can be accessed also via Bragg scattering P 
or through time of flight imaging, though indirectly since 
those techniques consist in measuring momentum distri- 
butions. Reconstruction of the full first order correla- 
tion function using absorption imaging, not yet obtained, 
would actually rely on measurements of the momentum 
distribution as a function of time of flight @. Studying 
the critical point of phase transitions and as a conse- 
quence fluctuations [9„-12j require precise spatial resolu- 
tion of correlation functions. In that trend, second or- 
der correlations have already been measured using Han- 
bury Brown- Twiss interferometry [13]. Accurate control 
of coherence properties will also be necessary for realizing 
quantum registers or quantum simulations with systems 
based on cold atoms trapped in optical lattices [14]. This 
control can consist in measuring phase correlations be- 
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tween lattice sites separated by Ax, i.e. the global coher- 
ence or s patially averaged f i rst orde r correlation function 
^(Ax) = ijj* (r) (r + Ax) /I?/; (r) p [l5|. This expression 
is true for a macroscopic wavefunction t/^, that may be 
approximated by ?/^(r) = \/iV^o(r) ~ ^d) for the 

example of a sample of N atoms in a ID optical lattice 
at sufficiently high lattice depths, with 0o(r) the macro- 
scopic profile [16| and w{x — nd) the Wannier function 
at site X = nd along the optical lattice axis x, d stand- 
ing for the lattice spacing. For a slowly varying profile 
in comparison to w and if |Ax| is an integral multiple 
of d one finds ^(Ax) = (j)^ (r) (/)o (r + Ax) because of the 
orthogonality of the Wannier functions. 

This article presents an alternative measurement of the 
global coherence function ^ in a Bose gas loaded in a ID 
optical lattice with a resolution of one lattice spacing. 
This study is based on superradiant scattering and em- 
ploys the mechanism of amplification of matter waves 
(MWA), that can be efficiently driven in our system as 
was shown in a previous experiment on mode compe- 
tition between superradiant scattering modes [T^- In- 
deed, following release of a ID optical lattice with period 
~ 0.5 /im formed with retroreflected laser beams, three 
different atomic momenta are mostly populated, the 
Bose gas at rest 0o(r) and two lightly populated waves 
a(j)o{r ± vAt) recoiling at =Fv velocity, with |a| << 1. 
The angle and frequency of a pump pulse for superradi- 
ance can then be chosen at the Bragg angle for optical 
diffraction (see Fig. [T]), so that the condition of am- 
plification of one of the wave is fulfilled [4]. Since the 
MWA gain directly relates to atomic coherence between 
the atomic reservoir and the recoiling seed, it is closely 
related to ^{vAt), where At is the time before the pump 
pulse for MWA is applied. Experimental analysis of this 
matter wave amplification process after release from the 
optical lattice shows the MWA gain evolves actually pe- 
riodically in time with a slow exponential damping due 
to lack of spatial coherence in the system. This result is 
used for characterizing coherence in the system while in 
the superfluide state. 

Our experiment (see Fig(T]) is performed with a cigar- 
shaped BEC of A^t - 2 X 10^ ^'^Rb atoms in the F = 
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Bragg mode 




time of flight perpendicularly to the x axis, from which 
we characterize the momentum distribution. Different 
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FIG. 1: (Color online) (top) Representation of the experi- 
ment. A cigar-shaped BEC is loaded along its long axis x 
in an optical lattice (OL) formed by two counter-propagating 
beams (Al ~ 852 nm), which creates a matter grating. Fol- 
lowing release of the optical lattice and a delay At, a 5 /is 
pump pulse (SR) is sent at ^ = 24° with respect to the x axis. 
This beam (Xs = 780 nm) can be superradiantly scattered 
into three modes, the Bragg mode symmetric of the incident 
beam with respect to the x axis and two end-fire modes along 
the X axis. (Bottom) Momentum representation of the Bose 
gas during the superradiance process. Only the main recoil 
atomic momenta are presented on the picture. Absorption of 
one pump photon by the condensate at rest (momentum 0) 
and subsequent emission at 24° into the Bragg mode leads 
to amplification of the matter wave (MWA) with momentum 
2hkL- The dotted circle stands for the energy conservation for 
atoms initially at rest during the Rayleigh scattering process. 



2,mF = 2 state, with Thomas- Fermi length 70 fim and 
width 7 /im [l8|, [l9|. The BEC, obtained after evapo- 
rative cooling in a QUIC trap, is loaded along its long 
axis in an optical lattice formed with a retro-reflected 
laser beam (A^;, = 27T/kL = 852 nm), focused on the 
BEC to a waist of 110 /im. A matter grating with period 
d = Al/2 forms because of the presence of the standing 
wave, whose optical lattice depth Vq which we are to tune 
from to 36 Er, Er = h^k\/2m being the recoil energy 
of a Rubidium atom. For an adiabatic loading to rela- 
tively high lattice depths, we use an exponential rising 
loading with 20 ms constant and 40 ms rise time, and 
then we keep the optical lattice switched on for 50 ms 
before sudden release of the combined optical and mag- 
netic traps. Following a short delay At, a light pulse, 
red-detuned from the D2 transition (A^ = 780 nm) by 
5 =1.3 GHz, of duration 5 /is and polarized along the 
z axis, is shined on the matter grating. The chosen an- 
gle « 24° between the incident light direction and the 
direction of the long axis of the condensate is given by 
the Bragg's law cosO = Xs/^l commonly used in X ray 
scattering. When this condition is satisfied, the pump 
pulse can be resonantly diffracted at 24° in a direction 
which is the symmetric of the incident beam one with 
respect to the long axis of the BEC as shown in Fig. [H 
Eventually, absorption imaging is performed after 30 ms 
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FIG. 2: (Color online) 30 ms time of flight absorption images. 
(A) Multiple matter-wave interference patterns after atoms 
were released from an optical lattice potential with a potential 
depth of 14.4 Er. (B) MWA for an initial loading at 14.4 Er 
depth and At = 0. (C) MWA at 36 Er lattice depth and 
At = 0. For sake of visibility, the upper scaling value in Fig. 
(C) is divided by 2. 



absorption images are shown in Fig. [21 for experiments 
performed with At = 0. In (A) is shown the multiple 
matter-wave interference patterns after atoms were re- 
leased from an optical lattice potential with Vq = 14.4 Er. 
In (B), when adding the pump pulse after release from the 
optical lattice, a strong enhancement of the -\-2JikL mode 
is observed because of MWA, while the and —2h\<iL or- 
ders are depleted. In (C), MWA for Vb = 36 Er appears 
still efficient, though the interference pattern is blurred 
because a complete loss of coherence occurs as the lat- 
tice depth is increased [f2lj. In principle, usual Rayleigh 
scattering superradiance (SRO) with light superradiant 
scattering into end-fire modes along the x axis could also 
be observed. Nevertheless, because of the large seeding 
of the +2fikL mode, SRO is completely negligible as a 
result of mode competition for lattice depths we consider 
in this article, i.e. Vq > 5Er [17]. 

For further explaining how to characterize coher- 
ence, interpretation of the data is carried out us- 
ing the quantum model of references IHI mod- 
ified to take into account coherence. In the su- 
perfluid regime, the one-particle wavefunction of the 
system at time t = immediately after lattice 
switch off is ip{r) = ^^^{x — nd)tljL{nd^y^ z) ~ 
V^L(r) (u;(0) +u;(-2kL)e2^^^^ +u;(2kL)e-2^^^^), where 
^l^L{'^) is the Fourier transform of the quasi- momentum 
distribution inside the optical lattice and w designates 
the Fourier transform of the Wannier function w{x). The 
matter wave field is limited to the first scattering orders, 
thus expressed as ^(r, t) = ^q(r, t)e*^'^Cq(t) summed 
over Co, 021^^ and c_2ki,, the annihilation mode oper- 
ators of atoms with corresponding momenta and pro- 
files ^q(r, t). Considering such profiles as approximately 
given by the free evolution of the initial wavefunction 
ip{r) during time t, we have (/)±2ki. (r, t) = (j)o{r ^ ^t^t) 
with (/)o(r, 0) = '0l(i*). Here, v = h2'kL/m is the velocity 
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at which the atoms recoil in mode 2ki,, related to the fre- 
quency of recoil uo = Ahk\/2m. The orthogonality rela- 
tion [cq, Cq/] ^ ^q^q/ Is Verified when the quasi- momentum 
spread 2A/c is smaller than 2/cl, A/c < /cl, which is valid 
if assuming only the lowest band is populated inside the 
optical lattice. The time evolution of the mode operators 
reads: 

C±2kz. = -i^C±2\^L ± ^ (GlC±2kz. + G2C^^2\^^ (1) 

In this expression, we have replaced cJcq by Nq the num- 
ber of atoms in the condensate at rest, that depends 
on time because of depletion during the superradiant 
process. As the number of atoms A^±2ki, of atoms in 
±2ki, modes are large, one replaces also the operators 
c±2kL by c-numbers of the form c±2kL = \/^V±2k^e*'^'^^^ 
with initial conditions c±2kL (0) = ^/A^SkT phase 
given by ip(t) = —uot during the free expansion 
time. Gi{t) = /dk|/3(k)|2|p(k,t)|2j[|k| - Ik^l] and 
G2(t) = /dk|/3(k)p|p(k,t)pe-^(^-^-+2J^-)-^(5[|k| - \ks\] 
are the gains for MWA into mode — 2ki,. They de- 
pend on the overlap integral p(k, t) = J (ir(/)o(r — 
vt, t)0o(r, t) exp[— z(k — k^ + 2kL) • r], centered at k^ — 
2kz., and on /3(k) = ^(^^-£^)i/2|k x z| the atomic 
coupling coefficient between pump and scattered light, ji 
and 1^0 referring, respectively, to the dipole moment of 
the transition involved and its Rabi frequency taken as a 
step- like pump pulse starting at time t = At. 

The possibility of retrieving the function g comes from 
the MWA dependence on the delay time At. Indeed, 
one can show the identity p(ks — 2ki,, At) = J dr(j)Q{r — 
vAt, At)^o(i*, At) = Q{vAt) which shows the gain for an 
exact transfer of 2JikL momentum, \l3{'ks)p0^s — 2kL)P, 
is not affected by the expansion of the profile (/)o(r. At) 
but only by the coherence Q{vAt). If one assumes 
^o(r, 0) = fo{y^ z)go{x)^ where /o is a gaussian func- 
tion independent on x or to the least varying very slowly 
with X and go depends only on x, the MWA gain of 
-\-2JikL mode, corresponding to an integration over the 
momentum distributions of the gains |/3(k5')p(k, At)p 
with k = ks, takes the form 

Gi{At) ^ A\P{ks)\^ J dxglix - vAt, x)gim{x - vAt, x) 

where A is a constant and gim{x — vAt^ x) = J dugi{u + 

X — vAt^u + x)e ^ , gi being given by gi{x — 
vAt^x) = g^ix — vAt)go{x) and >/2cr being the radius 
of the Bose gas inside the optical lattice, well approx- 
imated by the Thomas-Fermi radius in the superfiuide 

— tan^ 9q 

state. The Gaussian window function e that 
stems from the k^^ dispersion of the superradiantly emit- 
ted light, varies by less than 10% over 2a (our case 12 
lattice spacings d). Thus gim{x — vAt^x) provides a 
local measure of coherence at position x and taken in 
vAt when this window is larger than the local coher- 
ence length L. For the case of an ergodic system, the 



local coherence function is equivalent to the global co- 
herence, with consequence Gi{At) ex \Q{vAt)\'^. In re- 
ality, this approximation is not completely exact but 
even in case of sharp spatial variations of the state of 
the system with spatial variations of the local coherence, 
the function Gi still provides a good qualitative mea- 
sure of the global coherence. The variation of the am- 
plitude of the MWA gain with time At is then entirely 
contained in the Gi function, from which coherence be- 
tween positions separated by n lattice sites can be ob- 
tained when At = uTq with Tq = d/v = 7r/uj^4:0fis 
the time for a recoil over a distance equal to one lattice 
site d. G2 can be derived in the same manner, yielding 
G2{At) =A\(3{ks)\^J dxglix - vAt,x)gim{x,x ^ vAt). 
We can assume here gim{x^ x + vAt) ^ gim{x — vAt^ x) 
since the distances vAt that give a non zero gain are any- 
ways smaller than the coherence length, so that we define 
G = Gi ^ G2. Our model does not consider Doppler de- 
phasing during the very short superradiance pumping, as 
it appears in Eq. ([1]) with a unique recoiling frequency uj. 
Doppler dephasing occurring before pumping can be dis- 
carded too while deriving the expression for Gi as we are 
not interested in precision better than d in go{x). As can 
be shown rigorously, this is valid if the spatial spread 
induced by the velocity dispersion verifies VAkAt < 
which for the example of an exponential damping of co- 
herence with Ak 2/L only requires vAt < ttL. 

To verify this coherence dependence, we impinged a 
short MWA light pulse of 5 /is on the matter wave grat- 
ing after different delays At ranging from to 400 /is 
and for different lattice depths. From each time of ffight 
image, we compute the average momentum along x, 
(p) = J n{p^)p^dp^ ^ 2hkL{N2\,^ -7V_2kz.)/^^t. Such 
definition characterizes the average momentum change 
due to MWA even at high lattice depth, when differ- 
ent initial momenta are populated (cf Fig. [2] (C)). In 
FiglSl regular oscillations of {p) are observed as a func- 
tion of delay time for a 14.4 Er lattice depth. Damp- 
ing appears on a timescale of the order of 200 /is. To 
retrieve the damping time r and the oscillation period 
T, the data are fitted according to the heuristic formula 
A -h 5exp(-At/r) cos(27rAt/T + dcj)^), the initial phase 
5(j)Q being attributed to the finite duration of the pump 
pulse. We find T ^ Tq for every lattice depth. Such oscil- 
lations of the gain are directly linked to light diffraction 
by the matter grating, which is maximum only when the 
recoiling modes —2h'kL and 2h'kL are in phase. Accord- 
ing to Eq. ([T|), for small depletion of the condensate at 
rest and short pump pulses {Tp duration), the average 
momentum is directly proportional to G{At) and effec- 
tively oscillates with period Tq: 

(p) = 4hkLTp G{At) (1 + cos {2ujAt)) (2) 

We carried out a systematic measurement of the damp- 
ing time of oscillations versus the optical lattice depth 
(cf inset in Figl3]). Of the order of 200 /is in the super- 
fluid state, a fast drop of the damping time is observed 
between 15 and 20 Er, after which it is nearly zero. If 
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Delay time At (|as) 

FIG. 3: (Color online) (bottom) Average momentum (p) as 
a function of delay time At, calculated from time of flight 
images. The optical lattice depths are 14.4 Er (empty blue 
diamonds) and 36 Er (filled red circles). Error bars are ob- 
tained after averaging over several experimental data. Solid 
lines are empirical fits (see text). The black triangle points 
are a theoretical calculation of the global coherence function 
(reseated amplitude). (Inset, top) Dependence of the correla- 
tion length L versus the optical lattice depth obtained from 
the empirical fits of the data. Error bars indicate the 70 % 
confidence bound result from these temporal fits. 

the lattice depth is increased to 36 Er for example (see 
Fig El, almost no oscillation is observed, the damping 
time is indeed shorter than one period. From Eq. ([2]), 
we know this damping is directly related to the coherence 
of the system. We can infer a coherence length defined 
8iS L = vr to be larger than 5 d for Vq = 14.4 Er. If 
our interpretation, which is valid only in the superfluid 
state, could be extended to higher lattice depths, that 
would mean L is smaller than d for Vq > 17 Er. 

For comparing to theory, a coherence function propor- 
tional to |0('i;At)p is simulated for a 3D BEC of non 
interacting bosons loaded in a ID optical lattice in the 
superfluid phase, combining the models developed in ref. 
[23| and [2^ for the calculation of the system wavefunc- 



tion. Initial finite temperature of our BEC, T5 ~ 120 nK, 
is taken into account, so that limitation of the coherence 
comes from incoherent population of the different states, 
the magnetic trap levels being also considered for more 
precision. The computed coherence function is displayed 
in Fig. [3] for Vb = 14.4 Er. A relatively good agree- 
ment with the experimental results is observed in the 
superfluid phase. Smaller coherence at long distances is 
measured experimentally, probably due to the presence 
of interactions which should be considered in our sim- 
ulation. Our theory cannot account for the fast drop 
observed at high lattice depths either. This drop can 
be interpreted as a complete loss of coherence when the 
lattice depth is increased above a rather sharp thresh- 
old between superfluid and atomic squeezed state (21} . 
Large radial confinement at high lattice depths may also 
induce long-time scale redistribution of the atoms and 
consequently non adiabatic loading, heating and loss of 
coherence [25]. Further theoretical analysis will be nec- 
essary to determine the ultimate precision of our mea- 
surement of the correlation function and whether or not 
it is still valid above the observed threshold. 

In conclusion, MWA was used to characterize phase 
correlations of a Bose gas initially loaded in an optical 
lattice. A quantum model was developed to account for 
the main features observed in the experiment and con- 
firmed our ability to measure phase correlations with an 
accuracy of one lattice spacing in the superfluid phase. 
This method is relatively easy to implement. It may be 
used to characterize quantum phase transitions in well 
controlled systems [9|, for example when magnetic trap 
frequencies are relaxed to get less dense condensate and 
with 3D optical lattices. Superradiant light scattering 
rather than absorption imaging could be used as a direct 
way for measuring MWA gain. 
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